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Abstract 

We study stochastic billiards on general tables: a particle moves 
according to its constant velocity inside some domain T> C M d un- 
til it hits the boundary and bounces randomly inside according to 
some reflection law. We assume that the boundary of the domain is 
locally Lipschitz and almost everywhere continuously differentiable. 
The angle of the outgoing velocity with the inner normal vector has 
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a specified, absolutely continuous density. We construct the discrete 
time and the continuous time processes recording the sequence of hit- 
ting points on the boundary and the pair location/velocity. We mainly 
focus on the case of bounded domains. Then, we prove exponential 
ergodicity of these two Markov processes, we study their invariant 
distribution and their normal (Gaussian) fluctuations. Of particular 
interest is the case of the cosine reflection law: the stationary distri- 
butions for the two processes are uniform in this case, the discrete 
time chain is reversible though the continuous time process is quasi- 
reversible. Also in this case, we give a natural construction of a chord 
"picked at random" in T>, and we study the angle of intersection of 
the process with a {d — 1) -dimensional manifold contained in T>. 

Keywords: cosine law, Knudsen random walk, Knudsen regime, ran- 
dom chord, Bertrand paradox, kinetic equations, invariant measure, 
shake-and-bake algorithm 

AMS 2000 subject classifications: 60J25. Secondary: 37D50, 
58F15, 60J10. 

1 Introduction 

The purpose of this paper is to provide a rigorous mathematical treatment of 
a stochastic process that can be informally described as follows. A particle 
moves with constant speed inside some <i-dimensional domain. When it hits 
the domain boundary, it is reflected in some random direction, not depending 
on the incoming direction, and keeping the absolute value of its speed. One 
would like to understand basic stationary and dynamical properties of this 
process for "physically reasonable" domains and laws of reflection. 

Indeed, the current physics motivation for studying this process comes 
from the need to understand diffusive motion in porous media on a micro- 
scopic level. One wishes to explore the random motion of molecules inside 
pores of nanometer scale in order to deduce large-scale transport properties 
inside a porous grain. There are numerous applications of this basic problem. 
Recently, this has emerged as being of considerable importance e.g. in the 
design of synthetic zeolites for catalysis jTTj for which in a computer model a 
pore-selective molecular traffic control effect has been shown to enhance the 
effective reactivity of catalytic grains [6]. 

Inside a pore only very few molecules can reside so that the mean free 
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path between molecule-molecule collision is rather large. As a result intra- 
porous mobility is dominated by the interaction of molecules with the pore 
walls rather than among themselves. This is the so-called Knudsen regime 
in gas dynamics which gives rise to the model described above in which all 
molecule- molecule interactions are neglected [T8l [8] . As a further simplifying 
feature the kinetic energy of a particle is strictly conserved. Due to a com- 
plicated microscopic structure of the pore walls momentum is assumed to be 
transferred in a memoryless random fashion. 

Already Martin Knudsen proposed the cosine reflection law, in which case 
the reflected direction has a rotation-invariant cosine distribution around the 
surface normal. This is in some sense a physically natural choice (Feres and 
Yablonsky [TBI [Tl] ) and has attracted very considerable interest in recent 
years, see e.g. Coppens and collaborators [9l [TO], [2Tj , and also [26]. In partic- 
ular, it has been debated whether or not the self- diffusion coefficient (defined 
by the mean-square displacement of a particle) and the transport diffusion 
coefficient (defined via the density gradient between two open boundaries of 
a pore) coincide in the case of the Knudsen cosine law and how these diffusiv- 
ities depend on the structure of the pore walls. In two space dimensions the 
large scale motion of a particle inside a strip of constant width is analytically 
argued to be superdiffusive. 

The results reported in [9j [10l [2TJ [26] on diffusion in more complicated 
pore surfaces are mostly numerical simulation results for open pores where 
particles are injected and extracted. This numerical evidence, however, has 
proved to be subtle and often ambiguous, much depending on details of 
the definition of the quantities under consideration, which are frequently 
not provided. In order to clarify some of these open issues and for making 
further progress we provide here a mathematically rigorous treatment. In the 
present paper we focus on bounded domains which also prepares the ground 
for addressing the problems arising in future study of open domains. Pores in 
real solids such as zeolites have a rather complicated surface structure which 
is usually not known in great detail. Hence we keep our discussion as general 
as possible as far as the domain boundary and the reflection law is concerned. 
Nevertheless, particular emphasis is given to the previously studied Knudsen 
cosine reflection law. 

Next, we discuss related papers from the probabilistic literature. Moti- 
vated by an asymptotically minimax strategy for the so-called "princess and 
monster" game, stochastic billiards with the cosine law were considered by 
Lalley and Robbins [13 [2U] on a convex domain of the plane, and the invari- 
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ant measure for the position was found to be uniform. Later on, billiards 
have been used in Monte Carlo Markov chains algorithms by Borovkov [U [5] 
and Romeijn (25] (where the stochastic billiard with the cosine law of reflec- 
tion appears as "running shake-and-bake algorithm") to generate random 
vectors in a connected domain with smooth boundary or on the boundary 
itself. S. Evans [Uj gives a detailed account of the case of uniform reflection 
law and of a domain with C 1 boundary or a polygon. We emphasize that the 
present paper improves much on the geometry of the domain. 

From an analytical perspective, Bardos, Golse and Colonna pQ and Boatto 
and Golse [3] consider a gas of particles moving along straight lines between 
two infinite parallel horizontal planes, where the reflection angle is given by 
some dynamical system. They prove a diffusion limit in this model where 
there exists a natural angular cut-off. For some other models with determinis- 
tic behaviour inside the domain and stochastic behaviour on the boundaries, 
see Goldstein, Kipnis and Ianiro [15] and Caprino and Pulvirenti [TJ. 

We now briefly describe the content of this paper. It turns out to be 
useful to distinguish the random motion of the particle (below referred to 
as stochastic billiard) from the stochastic process defined by the sequence 
of hits of the particle at the boundary (below referred to as random walk). 
We consider domains which boundary is almost everywhere continuously dif- 
ferentiable and satisfies a Lipschitz condition. We first construct the walk 
and the billiard on such general domains. We then focus on bounded do- 
mains. For general reflection laws, we prove for both processes existence and 
uniqueness of an invariant probability measure, convergence to it, and a cen- 
tral limit theorem. A key observation is that these Markov processes satisfy 
the Doblin condition. 

The special case of the cosine law for the reflection, presents some interest- 
ing probabilistic, dynamic and geometric properties. The invariant measures 
are uniform, the random walk is reversible, as well as the stochastic billiard 
up to a sign change in the speed. The cosine law also yields a construction of 
a "chord picked at random" in the domain, which has universality properties 
independently of the domain. The mean chord length is equal, up to an uni- 
versal constant factor, to the volume-to-surface ratio of the domain. Also, 
for a smooth convex subdomain, our random chord induces a random chord 
on the subdomain, simply by its trace (conditionally on a non-empty inter- 
section): we prove that the law of this trace coincides with our construction 
of a random chord when directly performed on the subdomain. With a view 
on future applications to open systems we finally establish some statistical 
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property of the particle when it crosses a surface inside the domain. 

This paper is organized as follows. In Section [2] we define the model, 
comment on conditions for "good" behaviour and state our main results. 
Then, in Section [3] we establish some auxiliary results related mainly to 
geometric properties of domains with Lipschitz boundary. Our main results 
are proved in Section 0] (discrete time), Section (continuous time) and 
Section O (geometric properties of the random chord). 

2 Formal definitions and main results 

In Section 12.11 we describe the general conditions imposed on the domain T> 
and we introduce some notation related to this domain. In Section 12.21 the 
discrete time and continuous time processes are formally defined, and we 
discuss conditions under which the processes exist or have good behavior. 
The results are formulated in Section [231 

2.1 General notations and standing assumptions 

Let ||-|| be the Euclidean norm in IR d , and define B(x, e) = {y G M d : \\x— y\\ < 
e} to be the open e-neighborhood of x G M. d , and S ^ 1 = {y G : \\y\\ = 1} 
to be the unit sphere. Consider an open connected domain T> C M, d , and 
let &D be the boundary of T>, and T> be the closure of D (so that DU&D = T>). 
In this paper we will deal with three reference measures: 

• the ci-dimensional Lebesgue measure on T> or T>, 

• the Haar measure on the sphere § d_1 , 

• the (d — l)-dimensional Hausdorff measure in R d restricted to &D. 

To simplify notations we will use the same symbols dx, dv, dz, ... to denote all 
of them. We warn the reader, and we recall that ambiguity on the measures 
under consideration - should it arise - can be easily resolved by checking 
the space of integration. Similarly, when we write \A\, this corresponds to 
the c?-dimensional Lebesgue measure of A in the case A C T>, and to the 
(d — l)-dimensional Hausdorff measure of A in the case A C &D. 

Throughout this paper we suppose that dV is a (d— l)-dimensional surface 
satisfying the Lipschitz condition. This means that for any x G &D there 
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exist e x > 0, an affine isometry 3 X : IR d — > M. d , a function f x : M. d 1 — > R such 
that 

• / x satisfies Lipschitz condition, i.e., there exists a constant L x > such 
that — fx( z ')\ < As Ik — ^11 f° r an 2 ' (without restriction of 
generality we suppose that L x > 1); 

• 3 x x = 0, f x (0) = 0, and 

3 X (V n B(x, e x )) = {zE 13(0, e x ) : > f x (z^, . . . , 

Recall that, by Rademacher's theorem (cf. e.g. [121 theorem 3.1.6]), the Lip- 
schitz condition implies that the boundary dT> is a.e. differentiable. This, 
however, is not enough for our purposes, so we assume additionally that the 
boundary is a.e. continuously differentiable, and we denote by 1Z C dV the 
set of boundary points where &D is continuously differentiable. That is, we 
suppose that 1Z is open (with respect to the induced topology on dV), such 
that the (d — l)-dimensional Hausdorff measure of dT> \ 1Z is equal to zero, 
and dT> has locally C 1 parametrization in any point of 1Z. 

For all x G 1Z we can define a unique vector n(x) G with the following 
properties: 

(i) for all small enough £>0we have that x + en(x) G T>, 
(ii) 

inf (y — x) • n(x) = o(S) as 6 — > 

(i.e., n(x) is the normal vector which points inside the domain). Clearly, the 
map x —> n(x) is continuous for all x G TZ. 

Let x ■ y be the scalar product of x and y. For any u, v G we define 

by 

v) = arccos(w ■ t> ) 
the angle between u and v , and by 

S„ = {u> G S^ 1 : v ■ w > 0} 

the half-sphere looking in the direction v. If w G S> n (x), we use a simplified 
notation ip x (u) := </?(u, n(x)) G [0,7r/2). 
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Figure 1: On the definition of the quantities (p x (u), r x (u) and h x (u) 



Note that, by definition of n(x), if u G S n (x) then for all small enough t > 
it holds that x + tu G T>. Using this fact, for any pair (x, u) such that either 
x G 1Z, u G §„(£•)' or x &V, u & S ' -1 , define by (see Figured]) 

r x (u) = mi{t >0:x + tue dV} > 

the distance from x to &D along the direction u (we use the convention 
inf = +oo), and by 

h x (u) = x + ur x (u) G {dV, oo} 

the point on dT> seen from x in that direction (in fact, h x (u) = oo if and only 
if r x (u) = +oo). 

2.2 Construction of the walk and billiard 

In this section we define the main objects of interest, the continuous time 
random motion on T> described informally above, and its skeleton random 
walk, i.e., the discrete-time Markov chain corresponding to the sequence of 
points where the continuous time process hits &D. 
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Denote by e = (1, 0, . . . , 0) the first coordinate vector. Let r] = (r] n , n = 
1,2,...) be a sequence of i.i.d., 8 e -valued random vectors with absolutely 
continuous distribution on S e : there is a function 7 such that for any mea- 
surable B C § e it holds that P[r/ G B] = ^ B ^{u) du (where P stands for the 
probability). We suppose that 7 is supported on § e in the sense that 

inf 7(u) > (1) 

for any compact B C S e . 

To define the model, we need also a family (U x , x G dT>) of unitary linear 
operators in M. d with the property U x e = n(x) for all x G TZ, so that U x is a 
rotation of which sends e to n(x). Here, two cases need to be considered 
separately. 

• Case 1: the function 7 is rotationally symmetric around the first co- 
ordinate axis, i.e, 7(1*1) = 7(^2) whenever ui ■ e = u 2 ■ e. In this case 
we do not need any additional assumptions on the family (U x , x G dV) 
except measurability of x 1— > U x , and we denote TZ := 1Z. 

• Case 2: the function 7 is not symmetric. In this case we assume also 
that the map x — > U x is continuous for almost all x G &D, and we 
define TZ = {x G 7Z : the map x — > U x is continuous in x}. 

Fix an arbitrary x G 7?.; using the sequence (r] n ,n > 1), we define a 
discrete time Markov chain (£ n , n = 0, 1, 2, . . .) and a continuous time process 
(X t , t > 0), so-called stochastic billiard, in the following way. Below, (3 stands 
for a special state, that we allow for technical reasons, £ n or X t = & means 
that the corresponding process is undefined. 

(i) Put £0 = %o- Inductively, we define £ n+ i G {&D, 00, &} as follows: 

£n+i = < 6, if ^ = © or £n G <9D \ ft, 

[ 00, if £ n = 00. 

Since r] n has a density, one can show (see Lemma I3.2I below) that for 
any x G 7Z, the distribution of \\ x {U x rj) is absolutely continuous with 
respect to the (d — l)-dimensional Hausdorff measure on &D. Since 
almost all points of &D belong to TZ, almost surely it holds that either 
£ n G dV for all n, or there exists some n such that £ no = 00. 
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Figure 2: A stochastic billiard starting from the interior of V 
(ii) Again, we put X = x , and let us abbreviate r := and 



n = 1,2,... (formally, if ^ = 6 for some k, we put r k i = r k for all 
k' > k). Then, for t G [r n ,r n+ i) define 



Since a.s. the random walk never enters &D \ TZ, we obtain that, with 
probability 1, X t G T> is well defined for all t G [0, lim^oo r n ). 

Intuitively, the Markov chain £ can be described as follows: on each step, 
the particle chooses the direction of the jump according to the density 7 
"turned" in the direction of the normal vector (the one that looks inside the 
domain T>), and then instantly jumps to the point on the boundary of T> seen 
in this direction. As a process on the boundary of V, it is as close as possible 
to a random walk, and we will denote it like that. For some unbounded 
domains it may happen, of course, that such a point does not exist, in which 
case the particle jumps directly to the infinity. For the stochastic billiard X, 
when on the boundary, the particle uses the same procedure to choose the 
direction, but then it moves with speed 1 until it hits the boundary again (if 
it ever does). 

Let us note that we can define the stochastic billiard starting in any 
xi G V as well (in this case we have to specify also the initial direction), see 
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i=0 



x t = in + (t - r n )u^ n , 
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Figure [21 Of course, one may consider the random walk £ starting from any 
initial distribution on dT> (one has to require that the probability that the 
starting point belongs to 72. is 1); similarly, one may consider the stochastic 
billiard X starting from any initial distribution on T> x S d_1 (and, in fact, 
we will frequently do that in this paper). 

The next proposition shows that, supposing that the boundary of the 
domain satisfies the Lipschitz condition, the process X t is well defined for 
all t > 0. 

Proposition 2.1 With probability 1, lim n _ >0O r n = +oo. 

It would be interesting to consider more general laws of reflection: for 
example, one can suppose that the law of r\ is continuous, but is not supported 
on the whole § e ; and/or one can suppose that the law of rj has discrete and/or 
singular components. However, this would bring in some substantial technical 
difficulties. For instance, in the former case the domain may be such that 
some parts of the boundary are not seen from any of the boundary points 
(since we no longer have the rule "if the particle can go from x to y, then it 
can go from y to x"). In the latter case, there is an additional difficulty that, 
even if x G 72, it may happen that P[\\ x (U n ^ x )7)) G 72] < 1, and so one needs 
to define the set 72. differently (to ensure that with probability 1 the random 
walk never leaves the set 72), and show that, with that new definition, we 
still have \dV\TZ\ = 0. 

We end this subsection by commenting our assumption that the bound- 
ary satisfies Lipschitz condition. If one considers domain with the boundary 
which is not everywhere continuously differentiable (even if one supposes that 
it is not differentiable in only one point), then the Doblin condition ff23l) be- 
low need not be satisfied when the boundary is non-Lipschitz (for instance, 
it is not satisfied for the domain in Figure [3]). So, although it can be indeed 
interesting to study stochastic billiards on such domains, this task is much 
more difficult due to the lack of geometric ergodicity and related proper- 
ties (in [4j it was claimed that the results there are valid for domains with 
boundary consisting of a finite number of smooth surfaces; on the other hand, 
Theorem 1 of [I] is wrong for e.g. the domain of Figure [3]). Also, without 
Lipschitz condition it is not even clear if the process X t is well defined for 
all t > 0. If the density 7 is not symmetric, there are examples (cf. the 
one below) when the random walk £ n converges to a (singular) point on the 
boundary. It is, in our opinion, a challenging problem to find out if such 
examples exist for symmetric laws of reflection. 
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Figure 3: A stochastic billiard on a non-Lipschitz domain 



Example (Non-Lipschitz boundary). Consider the domain 

V = {(a, 6) G M 2 : < a < 1, -a 2 < b < a 2 } 

(see Figure E]). We denote by 

Gx = {(a, 6) G R 2 : < a < 1,6 = a 2 }, 
G 2 = {( a , 6) g M 2 : < a < 1,6= -a 2 }, 
G 3 = {( a , 6) G M 2 : a = 1,0 < 6 < 1}, 

the three smooth pieces of the boundary, so that &D = G\ U G 2 U G3. Let 
) G G 2 be a point on the lower part of the boundary, and 
abbreviate 9 X = arctan(2xi) (we use the same notation for x = (xi,x 2 ) G G± 
as well). Let Rp be the operator of anticlockwise rotation by angle (3 G [0, ir], 
and J be the symmetry operator defined by J(a, 6) = (—a, 6). The family of 
operators (U x ,x G dV) is then defined as follows: U x = R^ + g x J for x G G\, 
U x = Rs.-e x for x G G2, U x = J for x G G3. 
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Let g(x) be the first coordinate of i G I 2 . Define, for x G G\ (J G2, 

D x ((3)=g(K(Rpt(x)))-g(x), 

where t(x) is the tangent unit vector at x with positive first coordinate. 

After some elementary computations, one can obtain the following for 
small enough x: 

• D x ((3) < 1 for p < 8x1, 

• D X (J3) < *-f for/3e [&n;f + 9 X ], 

• D X (J3) <0for/?G + + 9 X ], 

• D X (J3) < - 2 -f for (3 > f + 0*. 

Consider the following density on § e : being /3 the angle between the 
element of S e and (—62), write (with some abuse of notation) j(/3) = A[/3 3 , 
(3 G (0, 71"). Since 

EG/fo+l) - <?(£n) I £n = X) = ^ D x (P)P 3 d{3, 

71 Jo 

one sees from the above above properties of D x (j3) that 

%(U)-3(Ulen^)<0 (2) 

for all small enough x. Define cr(?/o) = min{n : g(£ n ) > ?/o}- The equation (|2J) 
means that, for small enough t/ , the process g(£,nAa(y )) is a supermartingale. 
So, there exists = lim^oo g(CnAa(y ))j an d, by the Fatou lemma we have 
E(7oo < Eg(£ ). If £0 < Vo, this implies that cr(?/o) = 00 with positive 
probability, and then with a little more work one concludes that 

£ n — >■ a.s. as n — > 00. 

It seems clear that one should be able to construct examples of such kind 
where it holds also that lim^oo r n < 00. However in this particular example 
it is rather difficult to find out whether this happens, so we opted for not 
discussing it here. Anyway, existence of such examples is already a good 
reason to assume that the boundary of the domain satisfies the Lipschitz 
condition. □ 
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Figure 4: du = \\x — y\\ d 1 cos (3 dy 
2.3 Main results 

For the rest of this paper, we suppose that V is a bounded domain (that is, 
diam(P) < oo). Note that the Lipschitz condition on the boundary implies 
that in this case it holds that \&D\ < oo and that the number of connected 
components of &D is finite, see Lemma 13.11 below. 

For any x,y £ M, d , x ^ y, let us define the vector £ x>y = jnSr £ S* 2-1 ; 
note that £ x>y = —£y tX - We say that y £ T> is seen from x £ T> if there exists 
n £ § d_1 and £ > such that x + tu £ X> for all £ £ (0, t ) and x + £ M = V- 
Equivalently, the open segment (x, y) lies in V and x ^ y. Clearly, if y is 
seen from x then x is seen from y, and we write "x y" when this occurs. 

From the definition of the Markov chain £ n , we will show that the tran- 
sition kernel can be written 

P[£ n+ i £ A | i n = x] = / K(x, y) dy , 

J A 

where (see Figure Hj) 

K(x,y) = n — ,| rf , l{z,y £ U,x <-> n) (3) 

If — 2/11 

is the transition density. Formally, if x or y does not belong to 1Z, we put 
K(x,j/)=0. 

As mentioned above, we consider separately the important case of the 
cosine law of reflection 

7(u) = 7 d cosy3(e,w) = ld e ■ u, (4) 
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where the normalization constant 7^ is defined by 



jd — I cos if (e, u) du^j = ( / u-eduj 



making K( Xj 'J cL probability density on dV for all x E 1Z. The transition 
density then becomes 

K{x,y) = „ J|d _ 1 J l{x,y eTL,x ^ y} 



x-y\ 

.rn I f.T — 11) ■ nfv/n 

{x,y E K,X y} 



7d((y - g) - n (») (fo - y) • n (i/)) 1 j 



\\x — y\\ d+1 

When considering the specific case of the cosine law of reflection, we write 
£, X instead of £,X, and we call these processes Knudsen random walk 
(KRW) and Knudsen stochastic billiard (KSB). 

Define /to to be the "uniform" probability measure on &D: £io(A) = 
4^L. Since K(x,y) is symmetric, we immediately obtain that the KRW f 
is reversible, with the reversible (and thus invariant) measure /to- For other 
reflection laws it is usually not easy to find the exact form of the invariant 
measure (except for some particular cases, see [TTj) . but nevertheless we prove 
that such a measure exists and is unique, the random walk converges to it 
exponentially fast, and the Central Limit Theorem holds. 

Theorem 2.1 Suppose that diam(D) < 00. 

(i) There exists a unique probability measure fi on dT> which is invariant 
for the random walk^ n . Moreover, there exists a function ip '■ dT> — > R + 
such that fi(A) = f A ip(x) dx which satisfies 



i/)(x) = / ip(y)K(y,x)dy. (5) 

J&D 

Finally, the density ip can be chosen in such a way that inf,9x>?/> > 0. 

(ii) There exist positive constants Po,Pi (not depending on the initial dis- 
tribution of £0) such that 

||P[ae-]-/i|| v </5 e- ft ", (6) 

where \\ ■ || v is the total variation norm. 
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(Hi) We have the Central Limit Theorem: for any measurable A C dT> there 
exists a a, with a a > if < \A\ < \&D\, such that 

n 

i=i 

converges in distribution to a Normal(0, a\) random variable. 

In particular, for the KRW, §Df) and ffl) hold with ft = jx . 

We remark that the relation (jSj) for the Radon- Nyko dim density if> of the 
invariant measure fx is what is usually called the balance condition. 

Contrary to £ n , the process X t by itself is not Markovian, which, in princi- 
ple, makes it more difficult to analyze. To overcome that difficulty, we define 
another process Vt by 

Xi—r — Xt 
V t = lim— — -. 

ej.0 e 

Clearly, V tl = V 2 = U^Vn-i for all ti,t 2 G [r n _i,r n ). In words, V t is the 
cadlag version of the motion direction (or speed vector) of the stochastic 
billiard process at time t. As before, for the particular case of the cosine 
reflection law, we write V instead of V. 

Now, it is clear that the pair (X t , V t ) is a Markov process. Associated to 
this Markov process there is a whole family of martingales. Introduce the 
sigma-field $ t generated by X s , V s , s G [0, t], and C the transport operator 

Cf(x,v) = v ■ V x f. 

Let C 1,0 {T> x §> d ~ l ) be the set of functions on T> x § d_1 which are continuous 
and once continuously different iable in the first variable onPx S d_1 . 

Proposition 2.2 Let f G C lfi (V x S d_1 ) such that 

f(x, v)= f(x, U x u)"y(u)du Vx G 11, v ■ n(x) < 0. (8) 

Then, 

f(X t ,V t )- [ Cf(X s ,V s )ds (9) 
Jo 

is a martingale with respect to > 0). 
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From this, it follows that the generator of the Markov process (X t , V t ;t > 
0) is given by C on the set of smooth functions satisfying (jSJ). Note that the 
left-hand side of (jSJ) does not depend on v provided x &7Z,v ■ n(x) < 0. 

Next, we show that the stochastic billiard converges exponentially fast to 
equilibrium. 

Theorem 2.2 There exist a probability measure \ onT> x E> d ~ l and positive 
constants /3' , (3[ (not depending on the initial distribution of position and 
direction) such that 

\\P[X t e;V t e-}-x\\ v <Poe~ p,lt , (10) 

for all t > 0. 

We complement this result with a central limit theorem. 

Theorem 2.3 Let f : T> X — > R be measurable and bounded function. 
Suppose also that it is centered, that is, 

I f(x,v)x(dx,dv) = 0. (11) 
Then, as t — > oo ; the random variables 



'1/2 



f{X s ,V s )ds 



converge in distribution to a Normal (0, cr^) random variable, ah > 0. 

Remark. The variance ah may be equal to zero for some non-trivial /'s. 
Here is an example for KSB. Let G : T> — > R be a C l function, and 

f(x,v) = -v -VG(x). (12) 

Then, a"j = 0. The proof of this fact is placed in the end of the proof of 
Theorem 12.31 □ 



Next, we relate the invariant measure for the stochastic billiard to the one 
for the skeleton random walk. Recall that, by Theorem 12.11 the measure (x 
is absolutely continuous with respect to the Hausdorff measure dz on &D, 
and ijj(z) := ^(z). Let /jq and uq be the uniform measures on T> and S ' -1 
respectively, i.e., Ho(A) = j^j, u (B) = p^ry. 
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Theorem 2.4 The invariant measure \ is absolutely continuous with respect 
to ijlq <S> Vq, and is given by 



x(dx,dv) = ij}(z) z , . dxdv, z = h x (—v). 

COS(f z {V) 

In particular, the product measure /io <8> vo is invariant for KSB (X t , V t ) . 

For the rest of this section, we concentrate on the process with the cosine 
reflection law. 

By means of an usual procedure (shift the law of the KSB starting with 
the stationary measure /i <g> u by (— T), and take the limit as T — > oo) we 
can define a stationary version of the KSB for all t G R. 

Theorem 2.5 Let (X t ,V t ) t m be the stationary KSB, i.e., (X t ,V t ) is uni- 
formly distributed onVx S d_1 . Then, 

(X t ,V t ) teR '=(X_ t ,-V_ t ) teR . (13) 

This relation shows that KSB with uniform initial condition is not only 
stationary, but also reversible up to a sign change in the speed. Also, it yields 
a pathwise construction of the stationary version of the KSB for all i G R in 
the following way: first, choose a random point (xq, Vq) G T> x S d_1 according 
to the measure hq®vq. Construct a KSB (X^\ v} 1 ^ ) with the initial condition 
(X^\v^) = (xq,Vq), and an independent KSB (X^ 2 \v t ^) with the initial 
condition (X^ 2 \ V^) = (x , —v Q ). Then, the process (X t , V t )teR defined by 

f (Xi l \v t (1) ) fort>0 
(Xt, v t ) = I 

( (X (2) t ,-V^) fort<0 

is a stationary KSB. 

For any measurable A cV, define 

m(A) = E Ao / 1 l{io + U ioAl G A} dt 
Jo 

so xn(A) is the mean time that KSB starting from £o which is placed uniformly 
at random on dV, spends in A till the next hit of the boundary. In particular, 
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m := xn(T>) is the mean time interval between the consecutive hittings of the 
boundary. Use the abbreviation 



7r V2r(^±i)d 

m 1 r(f + l) 



ad! — x 1 



(for instance, k 2 = 7r, K3 = 4). Next, we calculate explicitly the quantity 
xn(A) (for the case of convex smooth domain, this formula can be obtained 
using the results of [1]). 

Theorem 2.6 Suppose that diam(P) < 00. We have 

m{A) = Kd W\' (14) 

In particular, m = K-dj^- 

The quantity m can be interpreted as "the mean chord length" for the 
domain T>, when the random chord constructed as follows: take a point on dT> 
uniformly at random, and draw a line from there using the cosine probability 
distribution. Formally: 

Definition 2.1 The random chord for a bounded domain T> satisfying our 
standing assumptions of Section \2.1\ is a pair of random variables (Hi,E 2 ) ; 
in dT>, with the joint density \dT>\~ 1 K(x, y). 

Note that, by symmetry of the above density, the construction defines an 
unoriented chord, as one expects from a random chord. 

Remark. Here is another construction of the random chord, which is equiv- 
alent to Definition 12.11 pick some (x, v ) in T> x § d_1 according to the den- 
sity p uS\\ {-v)\\ d x dv (C is the suitable normalizing constant), and put 
Hi := h x (v), H 2 := h x (—v). This is so because the KSB beginning with 
a uniform distribution on &D is a Palm version of the stationary KSB, and 
therefore the above fact follows from Theorem 4.1 of Chapter 8 of [28J and 
Theorem 12.41 □ 



With this definition, the quantity rh = E||Si — H 2 || is the mean chord 
length, and Theorem 12.61 has a nice geometric interpretation: the area (vol- 
ume, . . . ) of the domain is always proportional to the product of its perimeter 
(surface area, . . . ) and m: 

|*| = *W (15) 

K-d 



18 



This result also shows that one can obtain the quantity directly (i.e., 
without obtaining separately \T>\ and \&D\) by simulating the KRW £ n . 

Let us discuss some further properties of the above definition of random 
chord. The term "random chord" should remind the reader of the so-called 
Bertrand Paradox (see, for example, [16]), which can be briefly described as 
follows. A chord is chosen at random in a circle. The question is, what is 
the probability that the chord is longer than a side of the equilateral triangle 
inscribed in this circle? To answer this question, one considers three natural 
methods for choosing a chord at random. 

1. The "random endpoints" method: choose two points independently 
uniformly on the circumference and draw the chord joining these two 
points. Note that this is equivalent to choosing a uniform point on the 
circumference and drawing a chord from there with the angle (between 
the chord and the radius) chosen uniformly at random in (— 7r/2, 7r/2). 

2. The "random radius" method: choose uniformly a radius of the circle 
and a point on the radius and construct the chord whose midpoint is 
the chosen point. It is easily verified that this is equivalent to choosing 
a uniform point on the circumference and drawing a chord from there 
according to the cosine distribution, so this method is equivalent to our 
method of choosing a random chord on a circle. 

3. The "random midpoint" method: choose a point uniformly at random 
within the circle and construct a chord with the chosen point as its 
midpoint. 

By elementary computations, one finds that the above probability is equal to 
1/3, 1/2, 1/4 if the random chord is drawn from methods 1,2 and 3 respec- 
tively. It is then routinely argued, that the problem is not well formulated, 
since it is not clear which of the methods should be preferred, and the term 
"random chord" itself is not well defined. However, E. Jaynes [16] suggested 
that, among the three above methods, the method 2 is the "right" one, since 
it has certain invariance properties. (Recall that the method 2 is equivalent 
to our definition of random chord for the particular case of a circle.) In the 
next theorem we show that the random chord in the sense of Definition 12.11 
also has similar properties, thus leading us to suggest that perhaps Defini- 
tion 12.11 is the "right" definition of a random chord for general domains. 
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Figure 5: A random chord on V induced by a random chord on T> 

If T>' C T> is convex, a chord of T> which intersects T>' defines a unique 
chord on V by its intersection (see Figure [5]). Let us generate independent 
random chords (E 1 (i),E 2 (i)) oiV, i = 1,2,..., till the chord hits the do- 
main V, and then denote by (S'^S^) the intersection. We call (E'^E^) the 
induced chord on D'. Like in the "acceptance-rejection" algorithm (or "hit 
or miss") for random variable simulation, we easily check that (E'^E^) has 
the same law as the endpoints of [Hi, E 2 ] C\V' given that these sets intersect. 

Theorem 2.7 Let V (I'D be convex. Then, the chord induced on T>' by the 
random chord of T> is the random chord of T>' in the sense of Definition \2.1\ 

When the domain D' is not convex, the situation is more complicated, 
since a chord of D can intersect D' in several places. Denote by 



the chords obtained by intersecting D' with the chord (Ei,E 2 ) of D, with 
l G {0, 1, . . . , oo} their number. 

Theorem 2.8 Let D' C D be with the Lipschitz condition and an almost 
everywhere continuously differentiable boundary. For every Borel subset C 
of dD' x dD' , we have 




,k) > k < 1 



E 



\&D 



x \dV 



i 



L 



K(x, y)dx dy. 



(16) 



dD 



k=l 



In particular, 




\3V' 



(17) 



\dvy 
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Hence, the relation fll6p has the following interpretation: the expected value 
of the measure giving unit weight to each chord induced on V by the T>- 
random chord, is the product of the expected number of induced chords and 
the probability distribution for the (ordered) endpoints of the random chord 
of T>' . Another comment on this theorem is that the number i of induced 
chords is integrable, a property which does not seem easy to prove directly. 
The formula (fTTj) could remind the reader of the following fact concerning the 
so-called Poisson line process (see e.g. Section 8.4.2 of [27]): if V is compact 
and convex, then the number of lines hitting T> has Poisson distribution with 
the mean proportional to the perimeter of T>. 

Finally, we are interested in what happens at the moments when the KSB 
crosses a surface S C V. Suppose that S is a compact (but not necessarily 
connected) (d— l)-dimensional manifold (possibly with boundary) . Similarly 
to what we did for dT>, we assume also that S is locally Lipschitz and almost 
everywhere continuously differentiable. So, for cl. cl. X G S we can define a 
normal vector n(a;) -this time there are two possible choices- and consider a 
family of rotations (U x ,x G S) with the property U x e = n(x). We can make 
our construction so that the maps x i— > n(x) and x i— »• U x are continuous a.e. 
in S. 

Now consider the random surface crossing times f n of the KSB. We define 
this sequence of random variables r n , n = 1, 2, 3, ... in the following way: 

fi = inf{t >0:I ( G5}, 

and for n > 1 

f n+1 = inf{t > f n : X t G S}. 

It is not difficult to prove that for almost all the initial conditions (as well 
as e.g. when the starting location and direction are chosen independently 
and uniformly) the surface is intersected only finitely many times during any 
finite time interval: f n — > +oo asn-> 

Define 




u£v+ n , if y fn • n(xo > o, 

-Uf]V fn , if Vf n ■ n(X? n ) < 0, 



1 Indeed, first, analogously to Lemma 13. 2\ one can prove that a.s. all the points of 
crossing are regular points for S. Then, suppose that for some n the set {£„ + s(£ n +i — 
£n),s G [0,1]} fl S consists of infinitely many points. In this case, let yo G S be an 
accumulation point of this set; it is clear that ^„^„ +1 • n(yo) = 0. But then, as in the end 
of the proof of part (i) of Lemma l3.2l one can prove that this happens with probability 0. 
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i.e., w n is the relative direction in which the KSB crosses S for the nth time. 

By Theorems 12.21 and 12.41 we know that V t has asymptotically uniform 
(and independent of the location) distribution, so, if one observes the pro- 
cess when it is inside some set A C T> with positive Lebesgue measure, the 
empirical distribution of the direction should be asymptotically uniform as 
well. The next result shows that the situation is different when one considers 
the instances of intersection with a (d — l)-dimensional surface: in this case 
the law of intersection is actually the same as the law of reflection (this result 
is also connected to Theorem 12. 8j) . 

Theorem 2.9 For a. a. initial conditions and any measurable B C § e we 
have 



Remarks. 

• From the proof of Theorem 12.91 it can be seen that a stronger result 
holds: if we keep track also on the location of crossings, then in the long 
time limit, the location is uniformly distributed and is independent of 
the direction. 

• Even if one starts the KSB from the stationary measure, it is generally 
not true that (w n ,n > 1) are identically distributed, and, in particular, 
it is not true that w\ has the cosine law. Rather than constructing a 
counterexample formally, we refer the reader to Figure EJ it is easy to 
convince oneself that the law of w% should be close to uniform (and not 
cosine); more precisely, it converges to the uniform distribution as the 
number of vertical components of S converges to oo. 

3 Some auxiliary results and proof of the ex- 
istence of the stochastic billiard process 

First, we recall the definition of the a-dimensional Hausdorff measure. For 
any Borel set AcM'', define 



lim — y l{wi G B} = / 7dCos(^(e, u) du. 
n -*°° Jb 




(18) 




22 




Figure 6: With high probability, the (relative) direction of the first crossing 
of S is equal to the initial direction 



where tu a = 7r a / 2 /T(a/2 + 1), and the infimum is taken over all collections 
of open sets (Bj,j > 1) such that diam(Bj) < 5 for all j and A C {JJLi By 
Then, 

H a (A) = IimW?(A) 

<5|0 

defines the a-dimensional Hausdorff measure of the set A. 
We need to establish some properties of the domain T>. 

Lemma 3.1 Suppose that diam(P) < oo. Then we have \&D\ < oo and the 
number of connected components of dT> is finite. 

Proof. Let us show first that the area of a bounded surface satisfying the Lip- 
schitz condition is finite. There is a deep theory related to Lipschitz functions 
and Hausdorff measure that can be found in e.g. [12, 23J; nevertheless, since 
we only need this simple fact, we sketch its proof here. Consider the following 
elementary observation: if a function / satisfies the Lipschitz condition with 
the constant C then, for any bounded A C Z d_1 , 

diam{^ G Z d : (z (1 \ . . . , z^) G A, = f(z {1 \ . . . , z^)} 

< (C+ l) 1/2 diam(A). 

From the above definition of the Hausdorff measure it is elementary to obtain 
that the area of &D is locally finite; since &D is compact, our claim follows. 
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Let us prove that the number of connected components of &D is finite. 
Suppose it is not the case; then enumerate these components in some way and 
choose any point Xj from jth connected component, j = 1, 2, 3, . . .. Since T> 
is bounded, there exists x G V such that x is an accumulation point of the 
sequence (x i: i > 1). Since Xj G dV and dV is a closed set, we have x G dV. 
But then, in any arbitrarily small neighborhood of x one can find pieces of 
connected components different from that of x. This contradicts the Lipschitz 
condition, since we supposed that in B(x,e x ) there are no other connected 
components of &D. □ 

For any u G S ' -1 , a, b > 0, define 

£ u (a,b) = jrr G R d \ {0} : < a, tan tp(u, x/\\x\\) < &}, 

so <£u(a, 6) is a piece of a cone, looking in the direction u. 

The next lemma is a technical result that is needed, in particular, to show 
that the random variable £ is a.s. well-defined for all n. 

Lemma 3.2 Suppose that diam(X>) < oo. 

(i) Consider any xo G V and let B be a subset of positive measure ofS d ~ l 
such that, for some Eq > 0, we have Xq + ut G T> for all u G B, 
t G (0, £o)- Let ( be a B-valued random variable with the distribu- 
tion which is absolutely continuous with respect to the Haar measure 
on Then, the distribution of h xo (() ^ s absolutely continuous with 
respect to the (d—1)- dimensional H aus dor ff measure on dV. Moreover, 
with probability 1, ¥>h BO (o(4 BO (0,*o) < % / 2 - 

(ii) The conditional distribution o/£ n+ i given £ n = x\ G 71, is absolutely 
continuous with respect to the (d — 1)- dimensional Hausdorff measure 
on dV. 

Proof. First, we prove the part (i). We have to prove that the measure 
fi(A) = P[h XQ (Q G A] is absolutely continuous with respect to the (d — 1)- 
dimensional Hausdorff measure on &D. Consider any A C &D such that 
\A\ = 0. This means that, for any e > 0, there are sequences (zi G ~R d ,i > 1) 
and (rj > 0,i > 1) such that 

oo oo 

A c\J B( Zl ,r t ) and ^rf- 1 < e. 

i=i i=i 
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Without restricting the generality, one can suppose that the distance from x 
to any of B(zi,ri) is at least e /2. Let 

Bi = {u G § d_1 : there exists t > such that x + ut G B(zi, rj)} 

be the projection of B(zi, r^) — x on Clearly, there exists C £o > such 

that \Bi\ < Cegrf' 1 . If g(u) is the density of we can write 

P[h, (C)eA]< / g(u)du; 

on the other hand, 

oo oo oo 

I U < E E r ^ < 

i=l i=l i=l 

Since e is arbitrary, by the continuity property of the Lebesgue integral we 
obtain that ji(A) = 0. 

To complete the proof of part (i), observe that, with 7 the density of C, 

P[h X0 (()eA] = [ *,(v)i{h X0 (v)eA}dv 

= [ j(v)i{h xo (v) g Ann}dv 

= [ 7(4 „) ^%^ i{y g n,x - j,}dy. 
Ifo — y|r 

The second line is from the above absolute continuity, and the last line is 
obtained by the change of variables from v G S d_1 to y — h IO (f) G &D, 
which can be performed at regular points of the boundary provided that 

<fy(£y,x ) < 7T/2. 

Finally, it is an easy exercise to deduce the part (ii) from the part (i) 
(decompose § n (xi) into the union of Dj := {u G S n (xi) : r x 1 {u) e [V(i + 
l),l/j)}, j > 1, and D := {u G § n (a;i) : r x 1 { u ) — 1}> an d then use the 
part (i)). □ 

Define, for all x G &D, 

a x = sup{s > : there exists u G S d_1 such that x + C u (s, s) C V}. (19) 

Let us prove that, under the Lipschitz condition formulated in Section 12. 1[ 
the domain T> has the following properties: 
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Lemma 3.3 (i) For any compact A C M. d there exists a constant q(A) > 
such that a x > q{A) for all x G &D H A. 

(ii) Suppose that diam(P) < go. Then, for anyx G dV, there exists, e', 5 > 
andy x G TZ such that K(y, z) > 5 for all z G B(x, e)nTZ, y G B(y x , e'). 

(Hi) Suppose that x,y G TZ are such that K(x,y) > 0. Then there exist 
e, e', 5 > such that K(x', y') > 5 for all x' G B(x, e), y' G B(y, e'). 

Proof. Consider an arbitrary x G &D and fix oq < 1 in such a way that 
ao(l— Oq) -1 / 2 < (2L X ) _1 (recall that we assumed L x > 1). From the Lipschitz 
condition we obtain that for any y G dV H £>(£, ao^) it holds that y + 
^-^^((l - a o) £ x, L^ 1 ) C (TTli3(x, s^)). The collection (dV(~)B(x, a e x ), x G 
dT> fl A) is a covering of dT> H A; choose a finite subcovering to establish the 
property (i). 

Let us prove the property (ii). Suppose that diam(P) < oo and observe 
that one can obtain (see Figure [7]) 

V n B(x, e x ) n fl (y + £^ (ed) (oo, (2L,)- 1 )) ^ 0. 

yedT>nB(x,a 0£x ) 

This set being open, this implies that \B X \ > (here | • | stands for the 
(d — l)-dimensional Hausdorff measure), where 

B x = {z : \\z - x\\ = e x } (1 f| (y + € -i {e Joo,(2L x )- 1 )). 

ye&DnB(x,aos x ) 

Note that any point of B x is seen from any point of dV D B(x, a e x ), and for 
all y G B x it holds that 

^x{^x,y) > arctanL" 1 — arctan(2L 2 .) _1 . (20) 

Moreover, from Lemma 13.21 (i) it follows that one can find a point y x G TZ 
such that e x £ XiVx G B x , <p yx {£ y . X)X ) < n/2, and x <-> y x . 
For bi , b 2 > denote 

= £ M d : (z - y x - Wty^) ■ £ yxiX = 0,\\z- y x - < b }. 

Since y x G 7?. and (fy x (£ yx>x ) < tt/2, one can choose positive and small enough 
£>i, 6 2 , e' (also with the property bi < \\x — y x \\ — e x ) such that 
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• dVnB(y x ,e') C ft; 

• y <-» z for all z G G 6li62 , y e dV (1 B(y x , e'). 

Clearly, there is a small enough 62 < ^0 such that for the set 

H — {z : there exists t G [61, \\x — y x \\ — e x ] such that \\z — t£ yxjX \\ < b 2 } 

it holds that H d T> and {z £ H : \\x — z\\ = e x } C B x . Then, one 
can find positive and small enough s, e[ < e', ai < ao such that for any 
ye&DnBiy^ei) 

• (y + € iyx>x (oo, s)) n {z : \\z - x\\ = e x } C B x ; 

• (y + £w(°°> s )) n {z : (z - - My.,*) ■ £ Vx , x } c G 6o , 6l ; 

• ^(^oiea.) C (y + £ w (oo,s)). 

Then, by construction, we obtain that a?' «-> y for any x' G dV fl ai£ x )> 
y G <9X> fl B(y x ,e' 1 ). For part (ii) it remains only to prove that the above 
parameters can be chosen in such a way that K(x', y) is uniformly positive 
for regular points from the neighborhoods of x and y x . This can be obtained 
easily, using that (p y<c (&y x ,x) < tt/2 and (1201) . 

Proving part (iii) is very similar to the last part of the proof of part (ii), 
so we do not provide the details. □ 

Now, we are ready to prove that the stochastic billiard process X t is well 
defined for all t > 0. 

Proof of Proposition \2.1[ Fix any T > and let us prove that 

lim r n > T a.s. (21) 

n^oo 

Consider the compact set A = {x G lR d : \\x — Xq\\ < T} (recall that £0 = Xq)', 
by Lemma 13.31 (i), we have that a x > q(A) > for all x G &D fl A, where 
a x was defined in ([Til . Denote B u = {v G S d_1 : t&n<p(u,v) < q(A)}. Then, 
from (jTJ) it is elementary to obtain that 

P[||£n+i - U > a x \£ n = x}> inf P[r) G B u ) > (22) 

for any x G dT> fl A. This means that, as long as r n < T (and, consequently, 

£ n G dVnA), there is a uniformly positive probability that r n+1 — r n = 

£n|| > 9(^4) ■ This proves (l2"Tj) . and, since T is arbitrary, Proposition 12 . II holds. 

□ 
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4 Proofs for the random walk 



Here we prove the results concerning the discrete-time random walk £ n . 

Proof of Theorem \2.1\ First, our goal is to prove that the Doblin condition 
holds: there exist n , i > such that for all x, y G 1Z 

K no (x,y)>e. (23) 

To prove ( l23l . let us first suppose that x and y are in the same connected 
component of dT>. Now, from Lemma [3.31 it follows that for any x G &D we 
can find e(x), 5(x) > 0, s(x) G 71 such that 

inf K(y,z)>0. (24) 

z£llnB(s(x),6(x)) 

Note that (dV PI B(z, s(z)), z G dV) is a covering of dV, so let us choose 
a finite subcovering (dV n B(xi,e(xi)),i = 1,...,N). Abbreviate B i; j := 
dV n B(xi, e{xi)) n e(xj))> let 

j3 = min inf K(y,z), 

i=l,...,N y e_nnB{xi,e(x x )), 
z£llnB(s(x t ),6(xi)) 

Q\ = min I B t I , 

9 2 = min |9X>nB(s(a;i),<y(xi))| 

i=l,. ..,JV 

(note that all the quantities, defined above, are strictly positive). 

Since x,y are from the same connected component, there exists m < N 
and k(l), . . . , k(m) G {1, . . . , iV}, such that x G TZ n y e TZ (1 

B(xk( m )), and -Bfe(i),A:(i+i) 7^ for all z = 1, . . . , m — 1. Note also that one can 
suppose that m is the same for all x, y (and even that it does not depend 
on the connected component of &D to which belong x and y), since we 
formally can choose k(i) = k(i + 1) in the above sequence. So, abbreviating 
Dj = 1Z fl B(s(xk(j)), S(xk(j))), we can write 



K 2m (x,y) > j dz t I dw x I dz. 



dz m _i / dw m -i I dz 



Dm-l J -Bfc(m-l),fc(m) ^ Dn 
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K(x, z 1 )K(z 1 ,wi)K(w 1 , z 2 ) ... K{w m _2, Zm-l) 
x K(z m ^ 1 ,w m ^ 1 )K(w m ^ 1 ,z m )K(z m ,y) 

which proves (|23|) for the case when x, y are in the same connected component 
ofdV. 

Now, we sketch the proof of (1231) for the case when x, y are in different 
connected components of dV; we leave the details to the reader. Suppose 
that there are M connected components of &D, so that dV = U^Gj. Let 
dij = 1 if there exist G GifllZ, Zji G GjCiTZ such that K(z{j, Zji) > 0, and 

= otherwise. In words, A = (0^)^=1...^ is the matrix that describes 
the "graph" of the set of connected components of &D. It is clear that 
the matrix A is irreducible, let us prove that it is aperiodic. Let Gi be 
the connected component of dT> which is the boundary of the unique infinite 
connected component of M. d \V, we show that a ioio = 1. Let x be an exposed 
point of the convex hull of and let if be a supporting hyperplane of the 
convex hull of T> at xq, i.e., Hr\T> = {xo}. From the Lipschitz condition (see 
also the proof of Lemma [3.31 (i)) we know that for some u G S d_1 , b\, b 2 > 0, 
it holds that t u {b u b 2 ) C V. Let (see Figure E> 

H t = {tu + z : z G H is such that tu + az G T> for all a G [0, 1)}. 

Let 

e = inf{||x — y\\ : x G H,y G dV, \\y — x \\ = e x /2} > 0. 

Clearly, H t is a bounded set for all t G (0, eo) and H t R <9£> contains only 
the points of the connected component of Xq. Let Ci be a random variable 
with uniform distribution on (0, £0)5 and ( 2 a random variable with uniform 
distribution on S d_1 D (H — x ). Define 2/1 G such that 2/1 = du + a±( 2 
for some a\ > 0, and 2/2 £ -f^Ci sucn ^ na ^ ^ 2 = Ci M — a 2C2 for some a 2 > 0. 
By construction, we have that 2/1,2/2 are in the connected component of xo 
and 2/1 ^2/2- Similarly to the proof of Lemma [3T21 (i). one can show that a.s. 
2/1,2/2 G ft and ip(n(y 1 ),e yuy2 ) G [0,tt/2), y>(n(y 2 ), Zy 2 ,yi) e [0,vr/2). 

Thus, there exists m x such that A" 11 > 0. Consider any x E Gi, y £ 
Gj, there exist i = k ,ki, . . . , k mi _i, k mi = j such that = 1 for / = 

0, ... ,rai — 1. Then, roughly speaking, we go from x to z kokl in 2m steps, 

2 i.e., such that there exists an hyperplane H which intersects the convex hull of T> only 
at {x } 
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Figure 8: Finding two regular points which see each other 



then (recall Lemma [331 (hi)) jump from z fcofcl to z^ko, then go from z klko to 
Zkikz m ^ m steps etc. This shows that (123]) holds with no = 2m{m\ + 1). 

Now, having (1231) at hand, it is straightforward to observe that there 
exists a coupling of two versions £ and £' of the random walk with arbitrary 
starting conditions such that P[£ no = > e := e|9P|. This in turn means 
that one can construct the two processes £ and £' together with a random 
variable T in such a way that £f +fc = £^ +fc for all & > 0, and this random 

variable T (which is the exact coupling time for the two processes) has the 
property 

P[T > kn ] < (1 -e) k . 

Applying Theorem 5.1 of Chapter 4 from [28J, we obtain (jHj). 

The absolute continuity of fi with respect to the (d — l)-dimensional 
Hausdorff measure follows from Lemma [3.21 (ii). and the equation (jSJ) follows 
from the fact that if the distribution of £o is fi, then so is the distribution 
of £i. Note that i/j(x) = J gT> '4>(y)K no (y, x) dy for x G TZ, so from ( |23l) it 
follows that infgu ip > s. 

We end by proving (JTj). This is now completely standard, a direct conse- 
quence of (ii)-(iv), Theorem 17.0.1 in [22]. □ 
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5 Proofs for the continuous time stochastic 
billiard 

In this section we prove the results concerning the continuous-time stochastic 
billiard. 

Proof of Proposition ^. E For < u < t, we have 



E 



f(X t ,V t )-f(X u ,V u )- / Cf(X s ,V s )ds 



E 



(25) 



At, K n At) — /(^T„„iVu, K- n _i 



n>0 



TnAt 



Cf(X s ,V s )ds 



Tn-lVu 



5u 



^ E < r n , r„_! < t} (f(X TnM , K„At) - f(X Tn _ lVu , K- n _ lV «) 

Uli 



n>0 



TnAt 



Cf(X s ,V s )ds 

'T n _lVM 

X)e[i{« < r n < t}(f{X Tn ,VJ - f(X Tn , V Tn _ x 

n>0 

0, 



(26) 
(27) 
(28) 



where we have used Fubini's theorem in (|26|) . the fundamental theorem of 
calculus in (1271) . and the boundary conditions (JE]) in ( 1281) . This proves that ([9]) 
defines a martingale. □ 



Proof of Theorem \2.S[ To prove ( fTUl) . we use the approach of [28]. By 
compactness of the state space T> x § d_1 , there exists at least one invari- 
ant probability measure \ f° r the Markov process (X t ,Vt). Consider the 
process (X t , V t ) starting from an arbitrary initial position xq and initial di- 
rection v , and consider also the stationary version (X t , V t ) of this process 
(i.e., (X , Vq) ~ which implies that (X t , V t ) ~ x f° r an t). Our § oa l now 
is to construct an exact coupling of (X t ,V t ) and (X t ,V t ), i.e., we construct 
these processes together with a random variable T (the coupling time) on a 
same probability space in such a way that 

{Xf +a , V f+S ) = (X f+S , Vf +S ) 
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for all s > 0. Also, we show that T has exponential tail, which in its turn 
ensures ( TTU1) . 

Definition 5.1 We say that a pair of random variables (X, T), X G dV, T G 

K + is a- continuous on the set Ax B , A C dV, B C R+ ; i//or an?/ measurable 
A 1 dA,B 1 dB 

PLY eA u Te B x \ > a|Ai||Si|. 

For any x G dV, e > 0, define A^ = B(x, e) fl 7?.. In the sequel, we 
will write ^ for the position of the random walk that started in a, for 
the corresponding local times. Note that, if we consider the pair r"), its 
joint distribution is singular with respect to the Lebesgue measure (although 
typically the distribution of each component is absolutely continuous). The 
next lemma shows that in two steps the situation is already better. 

Fix an arbitrary x\ G 1Z with the following property: there exists yi G 1Z 
such that K(x\, y±) > and n(?/i) • £ yi)Xl < 1- 

Lemma 5.1 There exist a>i,e > 0, f 2 > f\ > such that, for any x G A S X1 , 
the pair (£f,Tf) is ct\-continuous on A£ x {f\,r<i). 

Proof. Choose small enough e, e' in such a way that K (x, y) > hi for all x G 
A^, y G A£ (this is possible by Lemma O (iii); note that A £ Xi U A e n C ^). 
Denote by 

= y G M d : y - 2/0 • n( m ) = 0} 
the tangent hyperplane at the point y±, and let 

vx = (n(yi) ■ ^^Jn^x) - £ 2/1)a . 1 . 

Note that 

• Vi + y% & H (in fact, t>! is the projection of the vector n^) — £ yitXl 
onto the hyperplane H — yi); 

• Vi ■ £ yuXl < (since vi ■ £ yuXl = (n(yx) ■ £ yi , Xl ) 2 - 1 is negative by the 
choice of yi). 

For notational convenience, for the rest of the proof we assume that d > 3; 
it is straightforward to adapt the proof for the case d — 2 (in fact, the proof 
in this case is much simpler, since the sets G, Gi, G^, G^ defined below would 
be simply one point sets). 



33 




34 



'J/1 1*1 



0. Let us 



Consider the (d — 2)-dimensional plane (see Figure [HD 

G = {z E R d : (2 - j/i) • Vi = 0} fl F. 

One can easily see that if y' — yi E G, then (y' — yi) ■ > 
choose an orthonormal basis t> 2 , • • • , i>d-i in G — yi, so that u 1, t> 2 , . . . , v^-i is 
an orthonormal basis in H — y\. Let V : 3D — > H be the projection operator 
defined by 

Vy = y-({y- yi) ■ n(j/i))n(yi), 
and denote, with -ff^y, 2) := ||x — + — z||, 

9 v (t) := R xl (V-\yi + ^),^i) = 2||P- I (y 1 + ut) -mil. 

As V^ 1 (yi + t>t) = y\ + vt + 0{t 2 ) when t — > 0, using also the fact that 
Vllxll = tAt, we have 



9' v (0) = 2 } yi + Vt - Xl }-v 

hi + vt-xtW 



t=0 



(29) 



(note that y x - x x = -£ yuXl \\yi -Xi||). 

Since t>i, . . . , t>d_i is a basis in H — y x , any s E H — yi can be uniquely 
represented as 



s = s (V + ... + s (<*-D^_ 1) 



so we regard . . . , s^ d ^ as the coordinates of s. Note that V 1 (s + y%) E 
&D and let us consider the function 



^ x , z (s) = R x (V- 1 (s + y 1 ),z). 
From fl29|) we obtain (recall that G is perpendicular to i Xl , yi ) 

^*xi,xi(s) 



9s« 



s=0 



s=0 



24i,yi • fi =: /i > 0, 
0, fori = 2,...,d-l. 



(30) 



Since *ff x ,z( s ) is continuously differentiable, for small enough e,e' we obtain 
that for x, z E A £ Xi and s + y\ E A^, 



3h • 
-,2h 



(31) 
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%M < * - / '• (32) 

Denote G x = {> £ G : \\z - yi\\ < e'/2}, G+ = G x + f^, Gl = G 1 - ^v x 
(again, see Figure [H]). Clearly, if e' is small enough, 

{s + tv x : a G Gi, t G [-e'/2, £72]} + ^ c P(A£). 

Let 

fi = sup sup *a,^(s), r 2 = inf inf ^^(s). 
Let us show that f x < r 2 . Indeed, 

MrO - M - ^) ^ t £ ' 

by (l3~Tj) . and for any s G G^ 



by (E2D, and the similar inequality holds for any s G G x . So, we obtain that 

f, f, \ 3ft _/ ft_/ h I _ h I 



Let 

M x ,M,r' 2 ) = {s G P(A^) — yi : ^ x , z (s) G (r[,r' 2 )}; 

observe that M XiZ (r^,r 2 ) = V' 1 {M x ,z{ r 'v r 2))- 

Note the following simple fact: if for all i£ [a\, 02] we have ci < g'(x) < 
C2 where < C\ < c 2 < 00, then for any / C [f(ai), /(02)] it holds that 
C2VI < I/" 1 ! 1 )! < c rVl- Using this for # s (t) = + v x t) (recall that, 

by ([31]), g'(t) G [3h/i,2h]) we obtain 



\M x ,M,r' 2 )\ > 



[ |^ 1 (K,r 2 ))|d S >|G 1 |(2/ i )- 1 K-r 2 ; 

JGi 



for all (r[, r' 2 ) C (fi, f 2 ). Clearly, for small enough e', there exists a constant d 
such that for any measurable A C "P^A^J we have |7 ?_1 (A)| > d\A\. This 
implies that there exists a' > such that for all x, z G A^ 

|M^(ri,r 2 )|>a'(r2-'"i)- (33) 



36 



Now, for arbitrary A 1 C A| , {r' 1: r' 2 ) C (f 1 ,r 2 ), an d a; G A| , using (1531) we 
can write 

PKfe A,rf e(r[,r' 2 )\ 



> 



A 

,2„y 



so (Cfi 7 "!) i s ai-continuous on A| x (fi,^) with «x = ^l 01 '- Lemma [5.11 is 
proved. □ 

We continue the proof of Theorem 12.21 Fix any b > 1 such that 

6 _1 f 2 — bfi > max{(6 — l)fi, (1 — 6 _1 )f 2 }. 

The following result is a logical continuation of Lemma 15.11 

Lemma 5.2 For any k > 1 t/iere exzsfo > stzc/j i/iai /or a// x e A| ; 

(^2fc? r 2fc) ^ s Oik- continuous on A| x (A;6fi, fc6~ 1 f 2 ) ■ 

Proof. We proceed by induction. The case = 1 follows from Lemma 15.11 
Let fc > 2 and suppose that (£%r k _u, T 2(k-i)) * s ttfc-i-continuous on A^ x ((k — 
l)6fi, (k — l)&f 2 ). Using Lemma \5. II and the induction assumption, write for 
any A\ C A £ X1 , (r[,r' 2 ) C (/cfofi, kb~ l f 2 ) 

^(fc-l)b _1 r 2 /■ /• 

> / / ak-idzdt / / ct\dz' di! 

J Ag 1 J(k-l)bfi J A 1 J(ri-t,r^-t)n(ri,r2) 

> Cajfe-iailAiK^-ri), (34) 
where 

min{(6 — l)fi, (1 — fr" 1 )^} 



min{f 2 — fi, (A; — l)(6 _1 f 2 ) — 6ri} 



One can obtain the last inequality in (together with the value of C given 
above) using the following fact: the convolution of functions C\l{x G [0, ai]} 
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and C 2 l{x G [0, a 2 }} (where a± < a 2 ) is a function supported on [0, a 4 + a 2 ], 
which grows linearly from to CiC 2 on the interval [0, a\], then is flat on the 
interval [ai, a 2 ], then decreases linearly from C\C 2 to on the interval [a 2 , ai + 
a 2 ]. So, if one considers any interval [c, c'] with c > and c' < ai+a 2 , then the 
convolution is uniformly positive (and one can easily calculate its minimum) 
on that interval. Note that in our case c = (b — l)fi, ai+a 2 — d — (l — b~ 1 )r 2 . 
Lemma 15.21 now follows from (j3"%|) . □ 

Now, let us choose fc in such a way that k (b~ 1 r 2 —br 1 ) > (n +l)diam(P), 
where n is from 0231) . Define 



^3 = kobr\ + nodiam(X>) , 
r 4 = k b~ l f 2 . 

It is important to observe that r 4 — r 3 > diam(P). 

Lemma 5.3 There exists a > suc/j t/iat (£^ 0+2fco , r ^ 0+2 A; ) ^-continuous 
on A| x (r 3 , r 4 ) /or any y £lZ. 

Proof. Consider any C A| and (r r , r") C (r 3 ,r 4 ). Let 3>^(t) be the 
distribution of r^ Q conditioned on = x}. Using (|23|) and Lemma [5T21 and 
the fact that the size of jumps of KRW cannot be larger than diam(P), write 

> / K n °(y, x) dx / P[C 2k0 e A,, rf fc0 G (r' - t, r" - t)] d*»(t) 

</A;s JO 



> eWJa^A^r" -r>), 

so Lemma |5~31 holds with d = e|A^ |afc . □ 

Now, we are ready to construct a coupling of two versions of the stochas- 
tic billiard process with different starting conditions. Let (X t , V t ) be the 
stochastic billiard with the initial condition (X ,V ) = (x ,v ) G V x 
and (X' t) V() has the initial condition (Xq, Vq) = (x' ,v' ). Let 

f = inf{t > : Xq + tv G <9£>}, u> = x + rt> 

be the time when (X t , Vt) first hits the boundary, and the hitting location, 
and let 

f = inf{t > : x' + tv' G dV}, w' = x' + f'v' 
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be the corresponding quantities for the process {X' v VI). 

By Lemma I5T3"1 the pairs of random variables (6n +2fc >^ + T n +2k ) and 
(Cno+2fe '^ / + r n '+2fc ) are k° tn a-continuous on A| x x (r 3 ,r 4 ). Abbreviate 
/i := r 4 — r 3 — diam(P) > 0. Since |f — f'| < diam(D), we obtain that 

\(r 3 + ?,n + f) n (r 3 + f',r 4 + f')l > ^ 

So, there exists a coupling of (f™ +2fco , f + C 0+2jfc J and (£™' +2fco , f' + r™' +2fco ) 
such that 

P[£i]>|A*>, (35) 

where 

-^1 = {£n +2fc = £n +2fc aIld ^ + T n +2k = ? + T n +2fc } 

(one can use e.g. the maximal coupling, cf. Section 4 of Chapter 1 of |28j). 
On the event E ll define T = f + T™ Q+2kQ . 

Suppose that the event E\ does not occur. Consider two cases: f + 

C +2ifco ^ f ' + r n '+2fc aIld f + T Z+2k >f ' + T no+2k - In tlle firSt CaSe > define 

™i = min{n > n + 2k : f + t™ > f + r™' +2jfco }, 

and m[ = n Q + 2k Q . In the second case, define 

m[ = mm{n > n + 2k :f' + T™'>f + r" +2fc(J }, 

^1 = no + 2fco (he., we identify the process with the smaller local time, and 
let it evolve a few more steps so that its local time becomes close to the other 
process' local time). Clearly, in both cases we have 

|(r + C)-(f' + r-;)|<diarn(P). 
Analogously, we obtain that one can couple (£™ l+no+2ko , f + r m 1+ no+2fc ) and 

(Cl+no+2fc ' f ' + T m[+n +2k ) in SUch a W that 

P[-&2 | E{\ > \ A £ Xi \h, (36) 

where 

^2 — {6mi+no+2fc = ^m' 1 +re +2fc aIld ^ + r mi+n +2fc = ^ + T m' 1 +n +2k } • 

On the event EfE 2 , define f = f + T™ l+TlQ+2ko . 
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Proceeding in this way, we construct two sequences of stopping times 
(m k , m' k ) k >i and a sequence of events 

rp / £W £Vj' pnrl -r -I- t w t 1 ' J- t w ' 1 

— lSm fc +no+2fe ~~ Sra' fc +n +2fco dI1U T ~r r m fc +no+2A:o — r ~r r m^+no+2fc J 

with the property (as in (155)) and (l3T)j) ) 

P[E fc |E 1 c ...^_ 1 ]>|A^ i |/ i ; (37) 
on the event E{... E c k _ x E k , we define f = f + T™ k+no+2ko . Since, 
C fc -r- fc _ i >2(n + 2/. )diam(P), 

and, by (1371) 

P[none of the events E%, . . . , E k occur] < (1 — | A £ Xl \h) k , 
we obtain that, for some (3[, /3 2 

P[f >t]<^e-^. (38) 

Up to time T, the processes X and X' can be explicitly constructed as in 
Section 12.21 Now, by definition we have that at time T the local times and 
the positions are the same for X and X', and so we have Xf = X'^. Since 
also Xf e &D, we can then construct the two realizations of KSB in such a 
way that Xf +s = X'^ +s for all s > (just use the same sequence of 77s from 

that moment on). This means that T is a coupling time, and, by Theorem 5.1 
of Chapter 4 from |28j, we obtain (flOl) from (138]) . This completes the proof 
of Theorem E21 □ 

Proof of Theorem \2.3[ The family of functions F t : T> x % d ~ x —>■ R, F t (x, v) = 
Jq E X)V /(X a , V,) (is is a Cauchy sequence with respect to the uniform norm 
|| • I |oo, in view of (ITT]) and (1TU1) . Hence it converges uniformly to a limit F, 
where 

poo 

F(x,v)= / E^/pC, V s ) ds. (39) 
./o 

Our interest for the limit is that 

M t = F(X t ,V t )- f f(X 8 ,V 8 )ds 
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is a martingale, a fact which follows from Q since F satisfies (JSJ) by definition. 
We can compute the bracket (M) of the martingale M (that is, (M) is 
a predictable process, uniquely defined by the following: (M) = 0, and 
M 2 — (M) t is a martingale) 

(M) t = *K < t}E x , v ([F(X Tn ,V Tn ) - F(X Tn ,V T -)] 2 | d T -) 

n>l 

Now, with N(t) > 1 defined by rjv(<)-i < t < r iv(t), we have 
— > a? as £ — > oo 



from Theorem 12. 1[ where 

o? = MgM^fMB e |o, oo). (40) 
E/iUCi — Coll 

By the central limit theorem for martingales [21 theorem 35.12], it follows 
that 

r 1/2 M t -> Normal(0, <rj) 
in law. Since F is bounded, we see that 

-l/2\ 



r 1/2 [ f(x s , v s )ds = -r^Mt + oit- 

Jo 



which implies the desired convergence and concludes the proof of Theo- 
rem [231 

Let us now prove that (|12p implies that a 2 = 0. For KSB, this function / 
is centered with respect to the uniform measure x which is invariant. In this 
case, we can compute the function F from fl39|) : 



F{x,v) = - / E X>V V S ■ VG{X s )ds 
Jo 



- lim E XtV / V. ■ VG(X s )ds 



G(x) - lim E x>v G{X t ) 
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= G(x) — \ G(x)x(dx 1 dv) 

which does not depend on v. Hence F(x,-) is constant on S n ( x ), and the 
definition (HOD shows that a) = 0. This can be understood in a different 
manner by computing 

M t = G(X t ) + [ V s ■ VG(X s )ds 



Therefore, 



G(x) — / G(x)x(dx 1 dv) 



[ f(X s ,V s )ds = G(X t )-G(x) 
Jo 



is bounded, which shows again that aj = 0. 

It is not clear to us if one can find a condition on / which is necessary 
and sufficient for a 2 f = 0. Notice however that if for a.e. x G &D, F(x, •) is 



constant on E> x , then a 2 f = 0. □ 



Proof of Theorem \2.4\ We identify the invariant measure x appearing in (fT0|) . 
By theorems 12.21 and I2.1[ we have for bounded continuous functions / : 
T> x — > R and starting from any initial condition, 

f f(x,v)dxdv = \imt~ 1 [ Ef(X s ,V s )ds 

JVxSd- 1 JO 

= hmEr- 1 ^ P f( X ^V s )ds 

n— >oo ' ' I 

k=l Jr k-1 

= (E A ||e! - eo||)~ 1 E A P ffa + sl^ds. (41) 

Jo 

The numerator in (l4"IT) is equal to 



<m 



/t(dz) / dv , y(U z 1 v) / dsf(z + st> , t> )l{z <-> 3 + st> } 



rfa; / ^||^-x|r d+1 ^(^)7(?7- 1 4,,)/(x,4,,)l{^a:} 



23 J923 <^ 
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with the change from polar coordinates (s, v) around z to Cartesian coor- 
dinates x. Then, by performing the change of variables from z G &D to 
v = £ z ,x G the Hausdorff measure integrator dz\\z — x\\~ d+1 becomes 

the Haar measure integrator dvj cosy? 2 (t>) with z = h x (— v), and the whole 
integral becomes 



dx 



v 



dz cos (fu x (-v){v) 



This is indeed the formula claimed in Theorem 12.41 



□ 



Proof of Theorem \2.5[ Since the process (X t , Vt)teR is stationary and Markov, 
it is enough to show that (X ,Vo, X t ,V t ) = (X t , — V t , Xo, Vq) for all t > 0, 
i.e., that 

E/(X , V , X t , V t ) = Ef(X t , -V t , X , -V ) (42) 
for all smooth test functions fonVx x D x S d_1 . Since 

Ef(X ,V ,X t ,V t ) = [f(X ,V ,X t ,V t )l{r n < t < r n+1 } 

n>0 

this equality (|42p will follow from the relation 



E 



f(X ,V ,X t ,V t )l{r n <t<r n+1 } 



E 



f(X t , -V t , X , -V )l{r n < t < r n+1 } 



(43) 



for n — 0, 1, The case n = being clear, we now consider the case n > 1. 

The left-hand side of ( ]4"3l is equal to 



dv dy 2 ... dy n+1 K(y 1 , y 2 ) . . . K(y n , y n+1 )f(x, v, y n + R£ yn , yn+1 , t yn , yn+1 ] 



where y\ = h x (v), where R is the function of x, y\, . . . , y n given by 

R = t- s 1>n - 1 1 7/i - x\\ , s kyn = \\y k+1 - y k \\ + . . . + \\y n - y n -i\\ , 

and A n C "DxS^" 1 x (dV) n is defined by \\yi— x\\+si )n <t< \\y\— a;||+si, n +i. 
Changing the Cartesian variable x into the polar coordinates s, w around y\, 
s G [0, 1 1 2/1 — yo\\], w G S d_1 (yo being defined as y = h x (— v)), we write 



dv dx 



dv\\yi-y \\ l{s < \\y x - y \\}dw ds 
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= -— dv -— K(y ,yi)l{s < yi - yo }os 

cosy2 yi (w) cos^ (tf) 

= dy Q dy x K (y ,y < ||yi - y ||}ds 

in terms of the Hausdorff measure on <92X Finally, the left-hand side of (1431) 
writes 

/ dj/o . . . dy n+ i dsK(y , y x ) . . . K(y n , y n+1 ) 

J R<\\Vn+l— 2/n||,S<||j/l— 2/0 1| 

with defined by + s + Si >n = £. Defining now 5 by S + r + s ljn = t, and 
changing variables to r = R, z = y n+ i, . . . z n+ i = y , the integral becomes 

/ dy ... dy n+1 dsK{z x , z ) . . . K{z n+1 , z n ) 

J S<\\z n+ i—Zn\\,r<\\zi—Zo\\ 

x / ( z n + S£ ZnjZn+1 , — £z n ,z n +u £yo,yn Zl ~ r ^yo,y\i ~^yo,yi) 

which is the right-hand side of (|43|) since K is a symmetric function. This 
ends the proof. □ 



6 Proofs of geometric properties of the ran- 
dom chord 

Proof of Theorem \2.(A To keep the paper self-contained, we provide the proof 
here, even though there are a lot of similarities with the proof of correspond- 
ing results in [31 ITU HI?]. 

Let us start the process from the uniform distribution on dT> (recall that it 
is the stationary distribution for KRW). Using fflUj) and the ergodic theorem 
one can write that 

-t 



pl = lim r 1 / 1{X S G A}ds 

\V\ t-*oo J Q 

lim r" 1 / l{X s e A}ds 







Jim (n- 1 Y^Ti - n-ij) V 1 [ ' HXs e A} ds 

i=i i=i ^ r »-i 
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= m' 1 m(A), 

so we already obtain that va(A) is proportional to \A\ (with a factor that 
does not depend on A): 

m(A) = ^\A\. (44) 

To calculate the value of the factor take any x G V and a small 5 > 
(so that £>(a;, 5) C V). For any u G S ' -1 denote by 

G u = {y G <9T> : there is s > such that y + su G £>(£, 5) and y ^ y + su} 

the projection of B(x, 5) in the (— w)-direction. For y G M. d let 

= |{s : y + su G B(x, 8)}\ 

be the length of the intersection of the ray drawn from y in the u-direction 
with B(x, 5). Then 

m(B(x,5)) = j^—r I dy ! j d cos (p y (u)H y>u du 

J dv J §n(y) 

/ du cos (p y (u)Hy iU dy. (45) 



Id 
\dV\ 



For almost all u G S d_1 there exists an unique y u G <9P that "sees x in 
direction it", i.e., {y u <-> x} and there exists i„ > such that y u + t u u = x. 
Define by 

G u = {z E M. d : (z — y u ) -£y U) x = 0, there is s > such that y + su G <B(x, 8)} 

the projection of i3(s, 5) onto the hyperplane that passes through y u and is 
perpendicular to £ VuiX - Then, provided that cos ip Vu (u) > 0, 

/ cos (p y (u)Hy tU dy = / (cos (p yu (u) + 0(5)) H VtU dy 

= I (l + 0(8))H z , u dz 

J G u 

= \B(x,5)\(l + 0(5)), 



so flU]) and ((45]) imply that 



m 

W\ 



Id 


\s d - 


-li 




dV\ 
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Figure 10: On the proof of Theorem 12.71 



and we finally obtain (fl4|) from (1441) . □ 

Proof of Theorem \2. 7[ First, we observe that the boundary of any convex 
domain is almost everywhere continuously differentiable. Indeed, since the 
boundary of a convex domain is locally a graph of a convex function, this 
follows from Theorems D and E of Section 44 of [23] . 

Consider a point y G dT>' where dV has a locally C 1 parametrization. As 
before, denote by 

G' u = {z G &D : there is s > such that z + su G dV' n B(y, 5),z <-> z + su} 

the projection of dV' n £>(?/, 5) onto <9D in the (— w)-direction. Let / = 
{[Si, S2] fl T>' 7^ 0} be the event that the random D-chord intersects T>' . For 
B C §n(y), write (cf. Figure [TU| and note that ^h 2) Hi = ^a',s') 

PpiGB(y,(5),^GS] 
= PfSieG^^s^GSl/] 

= PlSiGG^^H^GBlxPl/)- 1 
' c?t> / cosy^f) dx 



\dv\ P(j) 7 B 7 G: 



7d 



(|dl?'ni%,5)| +o(5 d - 1 )) / cos ^» (ft;. 
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This shows in fact that the couple (H^, ^w.h'J nas a density with respect to 
dydv which is proportional to cos f y (v )l{S n (j / )}, i.e., 



P[E' 1 eB(y,5),£ s>2 , s{ eB} 



Id 



\dV\ P(J) . 

which proves Theorem 12.71 
Remark. We can also identify the normalization in the last formula 



cos(p x (v) dxdv , (46) 
□ 



P(I) 



\dV'\ 
\&D\ 



(47) 



Proof of Theorem \2.8\ Let y G dT>' where &D' has a locally C 1 parametriza- 
tion, 5 > and a Borel set B C § n (y)- Similarly to fT46|) . we get 



3A;< t :5' eB(y,«), 



7/ =/ G 5 

- 2,fc>~l,fc 

Id 



l(dV'nB(y,8))xB 

Note that left-hand side can be written as 



cos <p x (v) dxdv. (48) 



E\^m, k eB(y,5),£ %A eB} 



k=l 



Then, for A\,A% G &D' such that A\ D A 2 = and x ^ y m V for all 
x G Ax,y G v4 2 , we get 



E 



K(x, y) dx dy. 



This is enough to ensure that, for C C dT>' x &D', we have 



k=l 



\&D\ 



K(x, y) dx dy. 



c 



which proves ([16]). Taking C = dV x <9£>' in (USD, we obtain (TTTj) . 
Proof of Theorem \2.9\ . For any v G S e and 5 > define 

= G M. d : there exist t G [—5, 5], a; G 5* such that z + tU x v = x} 



□ 
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Note that 

\S V (5)\ = 25(1 + o(5))\S\ cos <p(e,u). (49) 

Let us consider first the process in the stationary regime, i.e., we sup- 
pose that X , Vq (and hence X t , V t ) are independent and uniform. For 5 < 
va£x&&D,yes \\ x ~y\\ write (observe that the next probability does not depend 
on t) 

P [there exists n such that f n e [t,t + S],w n E B] 
= P[X t eSy t (5)] 

= -^Jp[X t eS u (S)\V t = u]du 

= 5(1 + o(5))\S\j^^ J^costp(e,u)du. 

This implies f|T8l) for the stationary case. For the general case, the for- 
mula ( fTHl) now easily follows from Theorem 12.21 □ 

Acknowledgements: We thank Frangois Delarue and Stefano Olla for stim- 
ulating discussions. 
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